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We investigate the radial behavior of galactic rotation curves by a Fourth Order Gravity adding 
also the Dark Matter component. The Fourth Order Gravity is a theory of Gravity described by 
Lagrangian generalizing the one of Hilbert Einstein containing a generic function of the Ricci scalar, 
the Ricci and Riemann tensor. A systematic analysis of rotation curves, in the Newtonian Limit of 
theory, induced by all galactic sub-structures of ordinary matter is shown. This analysis is presented 
for Fourth Order Gravity with and without Dark Matter. The outcomes are compared with respect 
to the classical outcomes of General Relativity. The gravitational potential of point-like mass is the 
usual potential corrected by two Yukawa terms. The rotation curve is higher or also lower than curve 
of General Relativity if in the Lagrangian the Ricci scalar square is dominant or not with respect 
to the contribution of the Ricci tensor square. The theoretical spatial behaviors of rotation curve 

<«_^ I are compared with the experimental data for the Milky Way and the galaxy NGC 3198. Although 

Cn , the Fourth Order Gravity gives more rotational contributions, in the limit of large distances the 

^ ' Keplerian behavior is ever present, and it is missing only if we add the Dark Matter component. 

Q^ ' However by modifying the theory of Gravity, consequently, also the spatial description of Dark 

^^ , Matter could undergo a modification and the free parameters of model can assume different values. 

After an analytical discussion of theoretical behaviors and the comparing with experimental evidence 
we can claim that any Fourth Order Gravity is not successful to explain the galactic rotation curves. 

Cn , In the last part of paper we analyze the gravitational potential induced by Lagrangian containing 

only powers of Ricci scalar. In this case we find an inconsistency in the boundary conditions in the 
fj ' passage from matter to the vacuum. 
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> ■ I. INTRODUCTION 



-Y~, ' Today the Universe appears spatially flat undergoing an accelerated expansion. There are many measurements 

If-, i: proving this pictures [ll-Q- According to the successful cosmological model [a-Qj there are two main ingredients in 
this scenario, namely Dark Matter (DM) and the cosmological constant A (Dark Energy). On the galactic scales, the 
evolution is driven by the usual Newtonian gravitational potential, but it needs hypothesizing the existence of DM 
to obtain a good experimental agreement. A good model for the galactic distribution of DM, in the framework of 
General Relativity (GR), is the Navarro- Frenk- White model (NFW model) [§|. 

However in recent years, the effort to give a physical explanation to the cosmic acceleration has attracted an amount 
of interest in so called Fourth Order Gravity (FOG), and particularly the /(i?)-Gravity, where / is a generic function 
of Ricci scalar R and the lagrangian of the theory is £ = f{R). These alternative models have been considered as a 
5— i ' viable mechanism to explain the cosmic acceleration. Other issues, of astrophysical nature, as the observed Pioneer 
. . , anomaly problem |10l . Illj can be framed into the same approach [l2l |. even if the anomaly has been shown to be 
attributable to well- understood classical mechanisms of radiation reaction [lj,[ij|- Apart the cosmological dynamics, 
a systematic analysis of such theories were performed at short scale and in the low energy limit |15l - [26| . 

Although FOG is investigated at theoretical level as alternatives to GR with Dark Energy and /or DM, an analytical 
approach to galactic rotation curves have been considered for /(i?)-Gravity [23, [231 ^^'^ ^^^ ^ i?"-Gravity [29l. [30{, 
with neQ, where the Hilbert-Einstein Lagrangian, C — R, is substituted by £ = R". In fact in these papers one 
does not consider a correction to scalar Ricci but one substitutes it by a power law for the Ricci scalar. In such way 
one find yet a Newtonian gravitational potential (ex r~^) with an added power law correction. 

By evaluating the rotation curve for a generic gravitational potential we want, in this paper, to use the potential 
induced by a /(A", F, Z)-Gravity, where for sake of simplicity we set X ^ R, Y = R^^Rap and Z = R^^'^^Rap^s- 
Then we generalize the Hilbert Einstein lagrangian by introducing a generic function depending not only on the Ricci 
scalar R but also on other two curvature invariants Y, Z where R^i, is the Ricci tensor and Rfi^ap is the Riemann 
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tensor. The gravitational potential is obtained when the Newtonian limit of theory with lagrangian C = f{X,Y,Z) 
is performed ,2^]. We are here interested to analyze the induced corrections when one considers the rotation curves. 
The corrections to potential can modify the global rotation curve and the DM model can undergo a modification. We 
report the radial behavior of rotation curve for two galaxies: Milky Way and NGC 3198. We analyze for any galactic 
substructure the rotation curve contribution in all frameworks: GR, GR with DM, FOG and FOG with DM. 

Since, today, there is not a successful model to explain the rotation curve without requiring the existence of DM, 
but there are many partially acceptable model, we conclude this work by comparing the outcomes of this paper with 
respect to the ones of i?"-Gravity. In this last topic we consider the analogy and the differences among the two 
theories by valuating also their metrics and Minkowskian limit. 

The plan of the paper is as follows. In Section [III we present the modified gravitational potential used. Section 
mil describes the rotation curve and the main properties induced in the potential. In the sections IIVI there is a little 
classification of principal galactic mass model and in section FVl we show the theoretical predictions of rotation curves 
and compare them with respect to the experimental data; section IVII resumes the analogies and differences between 
FOG and i?"-Gravity, while Section IVIII summarizes our principal conclusions. 



II. THE GRAVITATIONAL POTENTIAL BY f{X, Y, Z)-GRAVITY 

Let us start with a general class of FOG given by the action 



A = / d*xy/^ 



£ + X£„ 



= d x^f^ 



f{X, Y, Z) + XCr. 



(1) 



where / is an unspecified function of curvature invariants. The term Cm is the minimally coupled ordinary matter 
contribution. In the metric approach, the field equations are obtained by varying ([1]) with respect to g^j/. We get 



H = fxR + 2/yi?„^i?"'3 + 2fzRa.pysR'"'^^ - 2/ + n[3/x + fvR] + 2[(/y + 2/z)i?"'3].„^ ^ XT 



(2) 



-^, D = .^''^ and X = 8nG\ T.^ 



1 ^(v/^£„) 



is the energy-momentum 



where fx ^ ■^, fy = ■^, fz - jz^ ^ - ;- "-- - - ""- • -m- - ' ^^ ig^" 
tensor of matter and T is its trace. The second line of ([H is the trace of the first one. 

In the case of weak field and slow motion we consider the field equation in the so called Newtonian limit of theory. 
For our aim we can consider the metric tensor approximated as follows (for details, see [2ll [23. [3lll32| ) 



9ii^v 



l + 2$(i,x) 

-[l-2*(i,x)]^y 



(3) 



where $ and vp are the gravitational potentials and bij is the Kronecker delta. The set of coordinates^ adopted is 
x^ = (t, a;^, x^, a;'^) — (i,x). By introducing the quantities 



vn\ = 



m-i 



/x(0) 



3/xx(0)+2/i.(0)+2/z(0) 



/x(0) 



(4) 



/i'(0)+4/z(0) 



^ Here we use the convention c = 1. 

^ The Greek index runs from to 3; the Latin index runs from 1 to 3. 



we get three differential equations for the curvature invariant X and the gravitational potentials $, ^'^ 



(A-m2^)A$ 



(A - m.2'')R^, 



rn2 
2 



mi +2m2 A 



X = -vTia^A'p 



mi — m2 q2 
3mi'-i "ij 



m2 
2 



mi +2m2 ^ ^^ 



X = 



(5) 



I (A-mi2)X = TOi^A-p 



where A is the Laplacian in the flat space, Rij ~ A^ Sij + (^ — $),y is the ij'-componcnt of Ricci tensor and p is 
the matter density 25]. 

By choosing mi^ , ?Ti2^ > and introducing /11.2 = ■\/|toi.2^| the gravitational potentials in the case of point-like 
source {p = M S{x)) are given by 



'^Ak) = -^ 



1 + 1. e^'^il^l — - e^^^l^'l 



(6) 



1—1 g-Mllxl _ 2 g-/^2|x| 



If we have a generic matter source distribution p(x), $p/(x) becomes 



$(x) = -G I S 



, P(x') 



_ p-Mllx-x'l 



-P2|x-x'| 



(7) 



and an analogous relation is found for ^. The solution ([7]) has been obtained by using the superposition principle 
by starting from the solution ^. This approach is correct only in the Newtonian limit since a such limit correspond 
also to the linearized version of theory. The /(X, y, Z)-Gravity (like GR) is not linear, then we would have had to 
solve the field equations ^ for a given matter density. 

The parameters \Xi are the wave vectors but at same time in the fields theory are also the masses of propagation 
particles of field. In fact the equations ([5]) are the Newtonian limit (i.e. weak field limit and small velocity) of field 
equations ^ , while if we perform only the weak field limit of ([5]) we obtain the propagation of Ricci scalar (trace 
equation) and one of Ricci tensor (tensorial field equation) 



(n + m2^)i?^,y = source 
(n + TOi^)X — source 



(8) 



where D is the d'Alembert operator in the flat space. Then pi and p2 are the masses linked respectively to the 
propagation of Ricci scalar and Ricci tensor. 

If we choose the derivatives of / with respect to the curvature invariants satisfying the condition fxx (0) + Jy (0) + 
'^fziO) = wc find pi = p2 = p and the point-like solutions © become 



*p,(x) = *p,(x) 



GM 



-p|x| 



(9) 



then it is verified the condition gtt grr ^ — 1. A such condition is satisfied by the spherically symmetric metrics (for 
example Schwarzschild, Schwarzschild - dc Sitter, Einstein - de Sitter, Reissner - Nordostrom, etc). In fact in the 
paper [25| one has the condition A{^pi — ^pj) — ii pi — p2- Then we can affirm that only in GR the metric 
potentials $pi and 'i'pi are equals (or more generally their difference must be proportional to function |x|^^), while 
in FOG the Yukawa corrections must be equals. 



^ Throughout the paper we assume always fx{0) > 0, and therefore we may set /x{0) = 1 without loss of generality. 



III. ROTATION CURVES OF GALAXIES 

The motion of body embedded in the gravitational field is given by geodesic equation 



O-TT-jr = (10) 



ds2 "/^ ds ds 



where ds = ^Jgapdx^dx^ is the relativistic distance and F^o are the Christoffel symbols. In the Newtonian limit of 
theory we obtain from (1101) , formally, the classical structure of the motion equation 

d^~x- , , , , 

-^ - -V<I>(x) (11) 

but the gravitational potential is given by ([7]). A such potential is modified with respect to classical potential since 
we introduced other curvature invariants in the action ([T]). If we want to come back in the theory we have to set 
fxx = fy = fz = 0, then ^j.i , ^2 — >• 00 and the (O becomes the classical potential. 

The study of motion is very simple if we consider a particular symmetry of mass distribution p, otherwise the 
analytical solutions are not available. Our aim is to evaluate the corrections to the classical motion in the easiest 
situation: the circular motion. In this case wc do not consider the radial and vertical motion. The condition of 
stationary motion on the circular orbit is 



|x| 9|x| 



(12) 



where Vc is the velocity. 

The distribution of mass can be modeled simply by introducing two sets of coordinates: the spherical coordinates 
{r,6,(j)) and the cylindrical coordinates {R,9,z). An useful mathematical tool is the Gauss fiux theorem for Gravity: 
The gravitational flux through any closed surface is proportional to the enclosed mass. The law is expressed in terms 
of the gravitational field. The gravitational field g is defined so that the gravitational force experienced by a particle 
with mass m is Fgrav — mg. Since the Newtonian mechanics satisfies this theorem and, by thinking to a spherical 
system of mass distribution, we get, from (fT2t . the equation 



, ,2 GM(r) A-kG ,,.,,, 
vc{rf = ^ = / dyy'p{y) (13) 







where M{r) is the only mass enclosed in the sphere with radius r. The Green function of the /(X, Y, Z)-Gravity 
(^ |x — x'|~^), instead, does not satisfy the theorem [26|. In this case we must consider directly the gravitational 
potential fT]). Apart the mathematical difficulties incoming from the research of gravitational potential for a given 
mass distribution, the non-validity of Gauss theorem implies, for example, that a sphere can not be reduced to a 
point. In fact the gravitational potential generated by a ball (also with constant density) is depending also on the 
Fourier transform of ball [2g| . Only in the limit case where the radius of ball is small with respect to the distance we 
obtain the simple expression © . However in this paper we want to consider not the simple case of motion of body in 
the vacuum, but the more interesting case of motion in the matter. So we must leave any possibility of idealization 
and consider directly the calculation of the potential ([7]). 

Two last remarks on the ([7]) are needed. The two corrections have different algebraic sign, and in particular the 
Yukawa correction with pi implies a stronger gravitational force, while the second one (purely induced by Ricci and 
Riemann square) contributes with a repulsive force. By remembering that the motivations of extending the outcome 
of GR to new theories is supported by missing matter justifying the flat rotation curves of galaxies, the first correction 
is a nice candidate. A crucial point is given by the spatial range of correction. In fact the Yukawa corrections imply 
a massive propagation; then, more massive is the particle, shorter is the spatial range. In Fig. [T]we report the spatial 
behavior of gravitational potential ([5]) for arbitrary values interval of parameters pi and /i2. 

At last in Newtonian Mechanics the Gauss theorem gives us a spherically symmetric gravitational potential even if 
the spherically symmetric source is rotating. In GR as well as in FOG, however, the rotating spherically symmetric 
source generates an axially symmetric space-time (the well known Kerr metric) and only if the source is at rest 
one has the space-time with the same symmetry (the well known Schwarzschild metric). Then the galaxy being a 
rotating system will generate an axially symmetric space-time while we are using the solution ([7]). This aspect is not 




|x| 

FIG. 1: Plot of gravitational potential ((6|. fi2 = Cmi ^-nd fJ-i — .1 (dashed line), /ii = C/i2 and fi2 ~ -1 (dotted line) |25| ]. 
The behavior of GR is shown by the solid line. The dimensionless quantity ^ runs between -r- 10 with step 2. The dimension 
of /ii and /i2 is the inverse of length. 

contradictory because the solutions are calculated in the Newtonian limit (i.e. v^ <C 1) and under this assumption 
the Kerr metric collapses into Schwarzschild metric. In fact we have 
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(14) 



where r^ — 2GM is the Schwarzschild radius, S^ = r^ + rj^cos^O, H^ = r^ — r^r + rf , rj = L/M and L is the 
angular momentum along the z-axis. 

IV. MASS MODEL OF GALAXIES 

From the point of view of morphology, a galaxy can be modeled by considering at least two components: the bulge 
and the disk. Obviously the galaxy is a more complicated structure and there are others components, but for our aim 
this idealization is satisfactory. The bulge, generally, can be represented easily with cylindrical coordinates (but in a 
more crude idealization it is like a ball), while the disk has a radius bigger than the thickness. However we find that 
the rotation curve does not present the Keplerian behavior outside the matter, but the curve remain constant for any 
distance. Then we must formulate the existence of exotic matter that can justify the experimental observation. A 
simple discussion about the distribution of DM can be formulated by imposing the constant value of velocity in p^ 
for large distances. In fact we find 



Vc{r) 



constant 



PDAiir) ^ r 



(15) 



A matter distribution as (I15|) has a problem when we want to calculate the total mass. In fact if we have p ^ r^^ , 
the mass diverges. A such exotic behavior seems no-physical, but this outcome is only the consequence of constant 
rotation curve. In fact by increasing the distance also the mass must increase with the power law for any distance 
([T^. However, since the Gauss theorem holds in GR, the matter outside the sphere of integration does not contribute 
to the gravitational flux and then we do not have difference with respect to the ordinary matter. 

This same argumentation is not valid in f{X, Y, Z)-Gravity: the no viability of Gauss theorem implies that the 
range of integration of DM could cover all range and also the matter outside is considered. A cut off is needed now. 
In this paper then for completeness we consider that the galaxy is composed by three components: the bulge, the 
disk and an alone of DM. 



It should be noted that the spatial behaviors of DM (generally spherically symmetric) are made only a posteriori: 
the cornerstones of study of rotation curves are the GR and the distribution of ordinary matter. Only after this 
assumption the distribution of DM is such as to justify the gap between the theoretical prediction and the experimental 
observation. 

Before jumping to analysis of rotation curves we want to resume the principal spatial distributions of mass in the 
three galactic components. In literature there are many forms of density, but it is possible to resume them as follows. 

More realistic models are the ones with mass density depending also on the z-coordinate for bulge and disk. 
Particularly one can consider the following choice 
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(16) 
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:^a{ioM/is) r TT+^TIIF 



In this choice, suggested by Dehnen & Binney [3^, the bulge is described as a truncated power - law model (first line of 
(fT6| ) where Vb = /g dR' R' j^ dz'pbuige{R' ,z'). /3, 7, q and ^t are the parameters. While for the disk, one adopted 

iMo/pc^ and ^0 = 8.5 Kpc 






K+l 



a double exponential where the total mass is Md — Ztt^^ S0 e ^^ with E© = (48 ± 

[3J]. Finally in the case of DM the density profile is the NFW model [3,|3a] where g{x) = ln(l + x) 

^JJM ^'"^ '^'^ virial mass and virial radius and ^^ is a characteristic length. For the Milky Way £,'dmI^s = 
Leaving the axis-symmetry one can consider a more simple model: the spherical symmetry model, 
approach one has [36l - l38l | 
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where S,b, £,d, ^dm, M^, Md, Mdm, are the radii and the masses of bulge, disk and DM. k = 7.6695 and 77 = 22.665 
are dimensionless constants. The density profile of bulge considered is the well-known formula of de Vaucouleurs [39|. 
A more simple model, resuming the previous ones, can be 
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where r(x) is the Gamma function, < 7 < 3 is a free parameter and < a < 1 is the ratio of DM inside the 
sphere with radius £,dm with respect to the total DM. The radius £,dm and the mass Mjjm play conceptually the 
same role respectively of ^^'^ and M^^j. However, as before claimed, the hot point is the choice of DM model. Given 
all these models it is almost normal that there are many different estimates of DM. Therefore, the parameters of DM 
model may not be unique [40|. 



V. ROTATION CURVES BY f{X, Y, Z)-GRAVITY 

We are interesting to evaluate the circular velocity ([T^ adopting the mass models (US]) . So the potential ([7]) becomes 
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(19) 



where n is the distance on which we observe the rotation curve, ^ is the ehiptic function and the modulus of distance 
is given by 



A{R,R',z,z',e') = |x-x'| = ^{R + R'Y + (z ~ z'f - 4i?i?' cos' 



2 /)/ 



(20) 



The constant a is a scale factor defined by the substitution R ^r — >■ a r , a i? so all quantities are dimensionless. At 
last, by remembering r = Vi?^ + z^ the circular speed (|12p in the galactic plan is given by 



Vc{R) = \lR^HR,R,0) 



(21) 



In the Figs. [21 El we report the spatial behaviors of rotation curve induced by the bulge and disk component. The 
behavior for any component is compared in the framework of GR, FOG, GR + DM and FOG + DM. The values 
of free parameters of model are in the first line in Table U and refercing to Milky Way. The values of scale lengths 
/ii, fi2 are set at 10~^o~^, 10^ a^^. In both components we note for R >> ^b,^d the Keplerian behavior, while it is 
missing only when we consider also the DM component. The shape of the rotation curve is similar to ones obtained by 
varying the total mass and scale radius. For a given scale radius, the peak velocity varies proportionally to a square 
root of the mass. For a fixed total mass, the peak- velocity position moves inversely proportionally to the scale radius, 
or along a Keplerian line. 

As it is known in literature /(AT, 1", Z)-Gravity, and in particular /(A)-Gravity, mimics a partial contribution of 
DM. In fact the corrective term ex e~'^il''l/|x| contributes to enhance the attraction and thus the rotation curve must 
increase to balance the force. In the case of the other term, we have a correction ex — e~^^l''l/|x| that contributes, 
being repulsive, to decrease the velocity. However in both cases these terms are asymptotically null and f{X,Y,Z)- 
Gravity and GR must lead to the same result. Only with the addition of DM it is possible to raise the curve and have 
almost constant values. In the Fig. |4] we report the component of rotation curve induced by only auto-gravitating 
DM. 

In Fig. [5l we show the global behavior (experimentally expected) of rotation curve compared with respect to the 
bulge, disk and DM component for the /(A, Y, Z)-Gravity. While in Fig. |6] there is the global rotation curve in the 
framework of GR, FOG, GR -t- DM and FOG -I- DM. At last in Fig. [7] we replicate the outcome of Fig. Hbut we 
inserted the value /i2 = 5a^^. In this case the rotation curve induced by /(A', Y, Z)-Gravity allows lower values as 
previously we claimed. 
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FIG. 2: The rotation curve induced by bulge component (first line of HH)): GR (dashed line), GR + DM (dashed and dotted 
line), FOG (solid line) and FOG + DM (dotted line). 
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FIG. 3: The rotation curve induced by disk component (second line of H18[)): GR (dashed line), GR + DM (dashed and dotted 
line), FOG (solid hue) and FOG + DM (dotted line). 
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FIG. 4: The rotation curve induced by DM component (third line of fflE^ ): GR + DM (dashed and dotted line) and FOG + 
DM (dotted line). 



^ 

^ 



250 



200 



150 



100 



^^ 50 




R(Kpc) 



FIG. 5: Comparison between the rotation curves of galactic components: bulge (dashed line), disk (dotted line), DM (dotted 
and dashed line) and the global galactic rotation curve (solid line). All curves have been valuated in the framework of FOG + 
DM. 
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FIG. 6: The global rotation curve in the framework of GR (dashed line), GR + DM (dashed and dotted line), FOG (solid line) 
and FOG + DM (dotted line). The values of "masses" are ^i = 10"^ a"^ and /i2 = 10^ a"\ 
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FIG. 7: The global rotation curve in the framework of GR (dashed line), GR + DM (dashed and dotted line), FOG (solid line) 
and FOG + DM (dotted line). The values of "masses" are /ii = 10^^ a^^ and /i2 — 5a^^. 
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TABLE I: Parameters of models (O. The unity of mass is W^° Mq and a = 1 Kpc. 



Galaxy 



Mb Cb 7 Mi Cd Mdm io 



M Ot 



Milky Way 0.77 0.5 1.5 5.20 3.5 1.68 5.5 0.50 20 
NGC 3198 / / 2.60 3.5 0.84 5.5 0.53 20 



From the experimental point of view we used an updated rotation curve of Milky Way by integrating the existing 
data from the literature, and plot them in the same scale [41|. The data used are available in a digitized from 
the URL http://www.ioa.s. u-tokyo.ac.jp/~sofue/mw/rc2009/. The unified rotation curve shows clearly the three 
dominant components: bulge, disk, and flat rotation due to the DM p2 - |49| . These data, finally, have been updated 
further by [50| • The whole set of data are plotted in Fig. |8] and on them the theoretical rotation curve induced by 
/(X, Y, Z)-Gravity with DM has been superimposed. The values of best fit are shown in TableUwith /ii — 10~^ Kpc~^ 
and /X2 = 10^ Kpc" ^ 

The same mass model ([TS]) has been considered also for the galaxy NGC 3198. This galaxy has been chosen since 
the bulge is missing. Then we set Aff, — in the (|19p . In Fig. |9] we show the experimental data [5l| and the 
superposition of theoretical behavior. Also in this case we find a nice outcome for a new set of parameters shown in 
Table m while the values of /^i and ^2 are the same of Milky Way. 
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FIG. 8: Superposition of theoretical behaviors (GR (dashed hne), GR + DM (dashed and dotted line), FOG (sohd hne), 
FOG + DM (dotted line)) on the experimental data for Milky Way. The mass model used is shown in (|18p and the values of 
parameters are in Table|Il The values of "masses" are /ii = 10"^ Kpc~^ and ^2 ~ 10^ Kpc~^. 

By introducing the Gauss - Bonnet invariant defined by the relation Ggb = X^ — 4F + Z [54I we obtain in four 
dimensions the condition H'^^^ = Hf^^ — 4:H}^^ + H^^^ = 0. This condition on the field equations ([2]) enables us to 
consider only two curvature invariants [23], and by resolving the system ^ with respect to the quantities fxxiO) 
and /y(0) the general /(X, K, Z)-Gravity can be recast as the effective Lagrangian {Quadratic Lagrangian) in the 
Newtonian Limit 



c = /(x,y,z) = i?-3 



2^ll^ 



M2^ 



i?^ 






(22) 



The theory of Gravity represented by the Lagrangian (j22p is the more general theory considering all invariant cur- 
vatures, but we note a degeneracy. In fact we can have different /(X, y, Z)-Gravity describing however the same 
Newtonian Limit [2l|, [22, [2^. The solution of field equations or the experimentally detectable quantities, as the 
rotation curve, are parameterized only by the derivatives of /, then we can have different functions f{X, Y, Z) which 
admit the same physics. 

The initial aim, i.e. to extend the GR to a new class of theories, as we claimed in the introduction, is to justify the 
rotation curve without the DM component. From the previous outcomes, we see that even if the f{X, Y, Z)-Gravity, 
or better a /(X)-Gravity, admits a stronger attractive force, it is unable to realize our aim. Also in this framework 
we need Dark Matter. Obviously we need a smaller amount of DM on the middle distances, but for large distances 
we have the same problems of GR. 
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FIG. 9: Superposition of theoretical behaviors (GR (dashed line), GR + DM (dashed and dotted line), FOG (solid line), 
FOG + DM (dotted line)) on the experimental data for NGC 3198. The mass model used is shown in (|18p and the values of 
parameters are in Table IT] The values of "masses" are fii = 10~^ Kpc~^ and fi2 ~ 10^ Kpc~^. 

VI. 7?"-GRAVITY VS f{X, Y, Z)-GRAVITY 

The problem of DM seems to have been solved in literature, in the framework of /(X)-Gravity, by considering the 
Lagrangian C = i?" with n e Q [23, |30] . In these papers the gravitational potential for a point-like source can be 



$i?'.(r) 



GM 



1 



{r/rcY 



1 



(23) 



where re is a characteristic length and /3 is a dimensionless parameter. To recover the condition lim^-i-oo '^i?" ('') — 
one must have < ^ < 1. In the case /3 = the GR is found. 

We comment about the physical behavior of potential (P^ and we want to add some reflections considering the 
result of the rotation curve shown above. Before to analyze the mathematical properties of metric linked to potential 
(l^5|) . we want to show the different values of correction to the Newtonian potential. In Fig. [TU]we report the radial 
behavior of the corrections to l/r for the potentials ([6]) and ((23|) (to minimize the difference we considered only f{X)- 

- Gravity 



Gravity). From the plot we note a discrepancy between the two corrections. The correction by R 



6/ii 



acts over distances much smaller, while the correction induced by it!"-Gravity provides a potential nearly constant over 
large intervals and slowly goes to zero (~ r^~^). For this aspect the potential (P5|) does not need the DM component. 
Then with a procedure of fine tuning of re and j3 it was possible to justify the experimental rotation curve for a wide 
class of galaxies [30| when n = 3.5. This choice was possible because there must be a relationship (3 = f3{n) so that 
the potential (1231) w as compatible with respect to the field equations. These are the positive aspects of the potential 
([231) used in [s^i^l- 
We conclude this section by reviewing the fundamental weaknesses of i?"-Gravity. A) The potential ([^ presents 



an analogous behavior of potential ^ . In fact for r < 
like one induced by Ricci tensor square, while for r > 



one has 



(r/rc 



> then the correction is "repulsive" 



Tc one has an attractive correction. Now by remembering the 
reason of extension of GR, now we have unlike a repulsive contribution for r < re- If in f{X, Y, Z)-Gravity we can 
delete the Ricci tensor square contribution and we have only the /(X)-Gravity, in i?"-Gravity we must collapse only 
in GR. B) The potential (|23[) belongs to general class of solutions for i?"-Gravity classified by a perturbative method 
|53| . but the solutions are n-independent. Obviously the general solutions (it would be hard challenge to find them) 
are n-dependent, but at first order with respect to the perturbative parameter and in the vacuum ^ the field equations 
are identically vanishing. So we say that the presence of matter has been not considered and the choice of arbitrary 
constant has been evaluated only by matching _R"-Gravity with GR in the limit /3 — ?► 0. In fact by solving the field 
equations correctly in presence of matter (also with the point-like source) we would obtain solutions depending on 
the perturbative parameter and the technique is misplaced. 



This parameter is generally c ^, but the analysis is the same if we consider the dimensionless ratio rg/r. 
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FIG. 10: Comparison between the corrective terms induced by /(X)-Gravity (- 



solid line) and 7i!"-Gravity (- 
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dashed line), 
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/ii = 0.1, P — 0.8 and r^ = 5. The unities for r^ and /.ti are arbitrary. The dashed curve shows a very slow 



For these two aspects, but especially for the point B), i?"-Gravity does not admit the Newtonian limit if n ^ 1. 
The potential ((23)) does not follow a correct framework when extending the GR to the new theories we want to 
generalize the Newtonian potential. Generally all theories without Ricci scalar in the Lagrangian suffer from the same 
problem. For example also _R^-Gravity is in the same situation: is not possible to extend the solution in the matter 
[26| . Although we have solutions as 1/r with additional asymptotically flat terms, is not automatic the assertion that 
these solutions are the Newtonian limit of theory. 

Let us analyze now the mathematical properties of the metric trying to justify the difference of spatial behaviors 
in Fig. [TUl To simplify the calculation we choose a set of standard coordinates. The metric ([3]), from the expressions 
dH), becomes^ 



ds" 



1 4 

3" 3 



dt^ - 



1 + !1 ("l _ A^i^ + ^ e-^i^ - '2{fi2r+l) ^_^^A 



dr^-r^dn (24) 



where dVl = dO^ + siv? OdcjP' is the solid angle, while the element of distance linked to potential (|23)) can be written 
as follows 



ds^ = 



1 + 2$^"W 



dt^ 



l-2^fl^W 



dr^ 



'dil 



(25) 



where $^?(r) is the the potential (1^^ and ^^S'(r') is the other potential missing in the paper [30|- However in their 
analysis the knowledge of last potential is useless because its contribution in the geodesic motion is at fourth order. 
By following the paradigm of weak field limit at small velocity [2l|, l22| for the i?"-Gravity we find 



*i"l 



r) 



GM + K{I3) +Kx 13-1 GM 



1 

4r 



drr'^X{r) 



(26) 



where K{P) and Kx are constants depending, respectively, on the value of /3 and on the integral operation, while the 
Ricci scalar X could be an arbitrary function. In fact it needs some comment about the index n in the Ricci scalar. If 
n is a integer number, then the Ricci scalar can assume any value and can be also a generic space depending function. 
More attention is needed if n is a rational number. The field equations ^ take into account up to third derivatives 
with respect to the Ricci scalar, then we must ensure that the function f{X) and its derivatives are always well 
defined [5j| . In this case for n < 3 the solution Ricci flat {X = 0) or space depending and asymptotically vanishing 



The set of standard coordinates is defined by the condition to obtain the standard definition of the circumference with radius r. From 



the metric tensor ((3j we must impose the condition 1 — 2'^(r) r 



for the new radial coordinate. 
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are excluded. Only solutions with constant values are allowed, but the algebraic sign is crucial. A such behavior 
is expected any time we have the condition lirajf^o f{X) = constant [53|. Now in all these considerations we do 
not recover the condition linir-i-oo ^|;^(^) ~ 0: then we have a theory which do not provide us the Minkowskian 
limit. It is using the first perturbative contribution of a metric component (providing the flatness at infinity), while 
other contributions in the remaining metric components (negligible in the Newtonian limit) do not cover the same 
asymptotic limit. Then only for n > 3 we can have the flatness at infinity. 

Then we could say that for n > 3 the Minkowskian limit is recovered but a such perturbative approach can be 
performed only in the vacuum. The objection previously shown comes back. Up to third order (c~^ or Vg^) the 
geometrical side of field equation is identically null, but the matter side could not be null at first order (c~^ or rg). 

The class of it!"-Gravity are examples of theories where the weak field limit procedure does not generate automat- 
ically the Minkowskian limit. In fact only if we consider theories satisfying the condition limx-j-o /(^) = j53{ , 
their weak field limit is compatible with the request of asymptotically fiatness. Moreover /(X)-Gravity mimicking 
an additional source due to its scalar curvature [2l|, |5J] we would have a constant matter that pervades all space 
giving us a justification of more intense gravitational potential. In addition if limx-s-o f{X) = costant we do not 
have the Minkowskian limit, but we can interpret the apparent mass, only from the experimental point of view, as 
DM. These aspects, then, can be a mathematical motivation for different shape of point-like gravitational potential, 
but also source of further attention. 

VII. CONCLUSIONS 

In this paper we computed the study of galactic rotation curve when a FOG is considered. Among the several 
theories of fourth order we considered a generic function of Ricci scalar, Ricci and Riemann tensor. We started from 
the outcome of previous papers about the point-like solutions in the so-called weak field limit of theory and formulated 
the expression of rotation curve by inserting in the model the principal galactic components: the bulge, the disk and 
the DM component. In this limit we used the superposition principle for the potential since the Newtonian limit 
correspond also to the linearized version of theory. Since in FOG the Gauss theorem is not valid, the calculus of 
potential is performed directly by integrating in all space and obtaining in the potential also the information about 
the spatial shape of mass distribution. Apart the mathematical difficulties incoming from the research of gravitational 
potential the non-validity of Gauss theorem implies, for example, that a sphere can not be generally reduced to a 
point. 

Among two corrective contributions Yukawa-like to Newtonian potential seems that only one induced by a generic 
function of Ricci scalar is a good candidate; in fact this contribution has the same algebraic sign of Newtonian 
component and we can detect a more attractive force. A such situation can justify partially the bigger observed 
velocity, while the term induced by Ricci tensor (and by the Riemann tensor) acts lowering the velocity. 

The rotation curves have been evaluated by considering the bulge and the DM component spherically symmetric 
and the disk as a circular plane where the radius is larger than the thickness. Also in our case of FOG with ordinary 
matter we find that the rotation curve has the Keplerian behavior and only if we add the DM component we have 
a nice matching between theoretical and experimental data. However the hypothesis of existence of DM make two 
serious problems: since the matter distribution of DM is diverging when we consider the whole amount of mass is 
crucial the choice of cut-off inside the integral. Other hot point is the choice of the mass model of DM, or minimally 
the values of free parameters in the model. If we consider the GR as the theory of Gravity, we can have different 
values from those obtained if we choice the FOG. The spatial behaviors of DM (generally spherically symmetric) are 
made only a posteriori. 

The proposed models have been compared with the data of Milky Way and NGC 3198 obtaining the best fit for 
/ii = 10~^ Kpc~^ and /X2 = 10^ Kpc~^. In both cases we found a nice outcome for given values of free parameters of 
mass model. However the initial aim to extend the GR to a new class of theories, as we claimed in the introduction, has 
the hope to justify the rotation curve without the DM component. From the previous outcomes, while the f{X, Y, Z)- 
Gravity, or better a /(X)-Gravity, admits a more attractive force, on the other hand these modifications are not 
successful. Also in this framework we need Dark Matter. Obviously we need a smaller amount of DM on the middle 
distances, but for large distances we have the same problems of GR. 

We conclude the paper by comparing our outcome with the one of i?"-Gravity. In fact the problem of DM seemed 
to be solved by using a such theory. Particularly the constant rotation curves would be obtained without the DM. 
By finding a corrective term as power law in the radial distance one had a nice matching. But the model is not 
completely satisfactory. It is in our opinion that i?"-Gravity is misplaced but it contains some interesting and valid 
contributions. If the comparison with the experimental data encourages us to continue on this road, on the other 
hand, there are some mathematical pathologies that must be carefully evaluated. 

Furthermore, we think that a good weak field limit at small velocity, i.e. the Newtonian limit, of any theory of 
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Gravity must always satisfy the condition of Minkowskian flatness. In addition, a nice self-consistent theory, from 
a mathematical point of view, must be placed. It is obvious that can be considered also the weak field limit on the 
cosmological background and in this case we do not require the Minkowskian limit. 

Since, today, there is not a successful model to explain the rotation curve without requiring the existence of DM, but 
there are many partially acceptable models, we can say that, regardless of its detection, it remains an hard challenge 
to interpret the Dark Matter as a single geometric phenomenon. 
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